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Abstract 

The general Dirac equation in 1 + 1 dimensions with a potential with a completely general Lorentz structure is studied. 
Considering mixed vector-scalar-pseudoscalar square potentials, the states of relativistic fermions are investigated. 
This relativistic problem can be mapped into a effective Schrodinger equation for a square potential with repulsive 
and attractive delta-functions situated at the borders. An oscillatory transmission coefficient is found and resonant 
state energies are obtained. In a special case, the same bound energy spectrum for spinless particles is obtained, 
confirming the predictions of literature. We showed that existence of bound-state solutions are conditioned by the 
intensity of the pseudoscalar potential, which posses a critical value. 
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1. Introduction 

Since its formulation in 1928 [1], the Dirac equation has been widely investigated in various physical systems. 
Among the applications, we can highlight the Klein’s paradox [2-4], the description of electrons in graphene [5], the 
influence of the nuclear medium on the nucleons [6] and the relativistic hydrogen atom [7]. 

The pseudospin symmetry (PSS) is a topic of intense discussion, activities and recent progress in the last decades. 
PSS was introduced in nuclear physics for explain the degeneracies of orbitals in single particle spectra [8, 9]. In order 
to understand the origin of PSS, we need to take into account the motion of the nucleons in a relativistic mean field and 
thus consider the Dirac equation [10-12]. The case in which the mean field is composed by a vector {V,) and a scalar 
(Vj) potential, with 1, = Vt + Vs - 0 {Vi - -Vs) is usually pointed out as a necessary condition for occurrence of PSS 
in nuclei [10-12]. The study of symmetries in resonant states is certainly an interesting topic. The authors of Ref. [13] 
showed that the PSS in single particle resonant states in nuclei is exactly conserved, i.e., the pseudospin doublets with 
different quantum numbers ks -k + 1 have the same resonant state with energy and width Tres- That novel result 
were illustrated for single particle resonances in spherical square-well and Woods-Saxon potentials. Additionally, the 
Dirac equation exhibit spin symmetry (SS) when the vector and scalar potentials have the same magnitude and was 
used for explain the small spin-orbit splitting in hadrons [14]. Recently, a comprehensive review of the progress on the 
PSS and SS in various systems and potentials have been reported in [15], including extensions of the PSS study from 
stable to exotic nuclei, from non-confining to confining potentials, from local to non-local potentials, from central 
to tensor potentials, from bound to resonant states, from nucleon to anti-nucleon spectra, from nucleon to hyperon 
spectra, and from spherical to deformed nuclei. 

The four-dimensional Dirac equation with a mixture of spherically symmetric scalar, vector and tensor interactions 
can be reduced to the two-dimensional Dirac equation with a mixture of scalar, vector and pseudoscalar couplings 
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when the fermion is limited to move in just one direction {py = = 0) [16]. In this restricted motion the scalar and 

vector interactions preserve their Lorentz structures, while the tensor interaction becomes a pseudoscalar. This kind 
of dimensional reduction is very useful because the two-dimensional version of the Dirac equation can be thought 
as that one describing a fermion embedded in a four-dimensional space-time with either spin up or spin down [7]. 
Furthermore, the absence of angular momentum and spin-orbit interaction as well as the use of 2 x 2 matrices, instead 
of 4 X 4 matrices, allow us to explore the physical consequences of the negative-energy states in a mathematically 
simpler and more physically transparent way. Therefore, we can take advantage of the simplicity of the lowest 
dimensionality of the space-time. 

Considering mixed scalar-vector potentials the Dirac equation in (l-i-l) dimensions have been investigated for 
a sign potential [17] and smooth step potential [18]. In the context of mixed scalar-vector-pseudoscalar potentials 
(the most general Lorentz structure), the bound-state solutions for fermions and antifermions in (l-i-l) dimensions 
have been studied for harmonic oscillator potential [19], Poschl-Teller potential [20], Cornell potential [21, 22] 
and Coulomb potential [23]. In those works, the relation between spin and pseudospin symmetries using charge- 
conjugation and chiral transformation was illustrated. 

The square potentials, wells and barriers, are models widely used in low-dimensional systems such as the quantum 
dots [24], Dirac fermions in graphene [25], electrons in semiconductor heterostructures [26] and theoretical studies 
[27, 28]. Besides these applications, the well and barrier potentials are extremely examples used as toy-models in 
textbooks (for example, see [7] ), further increasing its importance in quantum mechanics. 

The main motivation of this paper is the approach of the Dirac equation (l-tl) dimensional in the framework 
of mixed vector-scalar-pseudoscalar square potentials. The scattering solutions furnish an oscillatory transmission 
coefficient, which does not present a total reflection. The presence of bound-state solutions are conditioned by the 
intensity of the pseudoscalar potential, which possess a critical value, compared to the mixed vector-scalar potential. 
Those interesting results are obtained as solutions of an effective Schrodinger equation for a square potential with 
repulsive and attractive delta-functions situated at the borders. The results can give us support for study PSS and SS 
in high dimensionality of the space-time. 


2. The Dirac equation in (l-l-l) D 


The time-independent Dirac equation for a fermion of rest mass m in the background of vector (V,), scalar (Vj) 
and pseudoscalar {Vp) potentials can written as (with units in which h - c - \ ) 


Hif/ - Elf/, 


H — CTip + cr^m + 


1 -H cr3 
2 




1 - 1 X 3 
2 


A -H 0-2 Vp 


( 1 ) 


where cri , cr 2 and 0-3 are the Pauli matrices and S = F, -H Vj, A = F, - Fj. 

The Dirac equation is covariant under x —> -x if Vp changes sign whereas F^ and F, remain the same. This is 
because the parity operator P - exp (is) Pocr^, where e is a constant phase and Pq changes x into -x, changes the sign 
of a-\ and <J2 but not of 0 - 3 . 

The charge-conjugation operation is accomplished by the transformation = o-ii/f* and the Dirac equation be¬ 
comes Hc^c - —Eipc, with 


He — cr\P + cr^m — 


I + cr 2 
2 


A- 


I — (Tt, 
2 


S + 0-2 Fp . 


( 2 ) 


One see that the charge-conjugation operation changes the sign of the energy and of the potentials Ff and Fp. In turn, 
this means that S turns into -A and A into -S. Therefore, to be invariant under charge conjugation, the Hamiltonian 
must contain only a scalar potential. 

The chiral operator for a Dirac spinor is the matrix = o"!. Under the discrete chiral transformation the spinor 
is transformed as - y^<lr and the transformed Hamiltonian Ply, - y^Hy^ is 


/ -H 0-3 I — CTj 

74 = 0-1/7 - 0-3OT -H —-— A -H —-— A + 0-2Fp. 


( 3 ) 


This means that the chiral transformation changes the sign of the mass and of the scalar and pseudoscalar potentials, 
thus turning S into A and vice versa. A chiral invariant Hamiltonian needs to have zero mass and F^ and Fp zero 
everywhere. 
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The equation (1) decomposes into two first-order equations for the upper, and the lower, i/f_, components of 
the spinor: 


dtj/- 

-i -h mtf/+ 

dx 


■ £(x)iA+ - iVpi//- = Eil/+ 


(4) 


.di^+ 

dx 


— + A{x)i//- + iVpif/+ — E[j/- 


(5) 


The components of four-current are given by and 7' = \f/^cr\il/. If we use the spinor \f/ in terms of its 

components the four-current is expressed by 7° = and 7' = 2Re(i/f* i/f_) which are conserved quantities 

for stationary states. Then, the Dirac spinor is normalized as dxJ^ = 1, so that i/r+ and are square integrable 
functions. It is clear from the pair of coupled first-order differential equations (4) and (5) that and have opposite 
parities if the Dirac equation is covariant under x —x. 

For A = 0 and E + -m, using the expression for obtained from (5) and inserting it into eq. (4) the Dirac 
equation for i/^+ becomes 


-A- = - 


I 


E + m 


T/ , 

HI - 


«•) 


dx^ 


, dVp 


4'^ 


- m^) l/r+ 


(7) 


Therefore, the solution of the relativistic problem is mapped into a Sturm-Liouville problem for the upper component 
of the Dirac spinor. As discussed in the ref. [20], we can take advantage of the discrete chiral transformation (y^) 
and we can obtain the solutions for S = 0 from the A = 0 case. This means that the chiral symmetry is invoked to 
obtain the equations obeyed by tf/+ and i/^_, for S = 0 and E + m. They are obtained from the previous ones by doing 
<-> if/-, m —> —m, S —> A and Vp —Vp. The solutions for A = 0 with E - —m and S = 0 with E - m, called 
isolated solutions [22, 29, 30], are obtained directly from the original first-order equations (4) and (5). For A = 0 and 
E = -m, we obtain 


d4f+ 

dx 


- Vpijf^ = 0 


( 8 ) 


whose solution is 


dip 

dx 


— -h Vpip- - —i (S -H 2m) if/+ 




(9) 

( 10 ) 


ip- = [N- - iNJ{x)] e-’'<^> 

where N+ and N- are normalization constants, and 

(11) 

v{x) = dyVpiy) 

(12) 

I{x) - J" dy p(j) -H 2m] 

(13) 


Note that this sort of isolated solution cannot describe scattering states and inasmuch as ij/+ and are normalizable 
functions, the possible isolated solution implies that Vp + 0, therefore the presence of a pseudoscalar potential is sine 
qua non for provide isolated solutions [22]. 
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3. The Square Potentials 

Lets us consider 


I.(x) = C^gix), 


(14) 


A = 0, 


(15) 


Vp(x) = Cpgix) 


(16) 


gix) = - [sgn(x -a)- sgn(x + a)] 


(17) 


where sgn(x) = x/|x| (x 0) is the sign function, and Cp are constants with dimensions of energy. Due to the 
chiral symmetry we can focus the discussion on the L case (A = 0). The results for the case A 0, X = 0 and Vp, still 
given by (16), can be easily obtained by just changing the signs of m and Cp in the relevant expressions. 

For the potentials (14), (15) and (16) we have not found a normalizable isolated solution with E - -m, i.e., the 


absence of isolated solution is because Vp{x) - 0 for |x| > a and therefore is a constant. For the case E + -m the 
Eq. (7) takes the form 



where 



5{x) is the Dirac delta function and the effective energy particle is given by j2m, with - E^- m^. We note that the 


effective potential distinguishes particles and anti-particles if + 0, and so we can not expect a symmetric spectrum 
with respect to E. 

Let us introduce the parameter r = Veffd-^l < a), where 


Veff(|x| <a)^Cl- C^(E + m) . 


( 20 ) 


The parameter r characterizes three different profiles for the effective potential as illustrated in figure (1). If r < 0, 
the effective potential consists in a finite square well potential at the region -a < x < +a with attractive and repulsive 
delta function situated at x = -a and x = +a, respectively. If r = 0, the effective potential consists in a double delta 
function potential with attractive and repulsive delta function situated at x = -a and x = +a, respectively. Finally, if 
r > 0, the effective potential consists in a finite square barrier potential at the region -a < x < +a with attractive and 
repulsive delta function situated at x = -a and x = +a, respectively. 

3.1. Scattering States 

We focused our attention to the scattering states solutions that describes a fermion moving from left to right. In 
this way, ^(x ^ -oo) describes an incident wave moving to the right and a reflected wave moving to the left, and 
ifi{x -oo) describes a transmitted wave moving to the right or an evanescent wave. The upper components for 
scattering states are written as 



( 21 ) 


where 


j] — — r 

The group velocity of the waves described above is given by 


( 22 ) 


dE k 



( 23 ) 







Figure 1: Profiles for effective potential Veff(x) for r > 0 (pointed line), r = 0 (thin line) and r < 0 (dashed line). 


where the double signal is related to propagation direction. For both range the probability current densities are given 
by 


= --(|AJ2_|A^|2) for E^-m (24) 

E + m 

and 

2k 

j\x +oo) = +- \D±f' for E ^ -m (25) 

E + m 

Note that /'(-cm) = Jmc-Jref and 7'(+oo) = Jtrans, where Ji„c, Jref and Jircms are nonnegative quantities characterizing 
the incident, reflected and transmitted waves, respectively. If £ + m > 0, then (A_e“'*') will describe the 

incident (reflected) wave, and D - 0. On the other hand, if £ + m < 0, then A_e (A+e^'^^) will describe the 
incident (reflected) wave, and D+ - 0. 

To determinate the reflection and transmission coefficients, we use the probability current densities given by (24) 
and (25). The ^-independent probability current allow us to define the reflection and transmission coefficients as 

R± = with R^ + T^^l (26) 

where the quantities 

= — and — (27) 

A± A± 

are called reflection and transmission amplitudes, respectively. 

We demand that i/r+ be continuous at x - ±a, that is 

lmtj,{x)\lztt>0. (28) 


Moreover, the effect due to the delta function potential on dtj/^ldx in the neighborhood of x — ±a can be evaluated by 
integrating (18) from +a - cr to +« H- cr and taking the limit cr ^ 0. Thereby, we obtain 


lim 

cr-^O 


d>j/+{x) 

dx 


\x~±a+<T _ 
'X=±a~<T~ 


Cp 

±—‘A+(-^ = ±a)- 


( 29 ) 
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With 4>+(-X) given by (21), conditions (28) and (29) imply that 


+ A_e 


Aka 


(30) 


= B+e""' + 


C„ 


- iT][B+e'''‘‘ - 8^6-'''“) = (B+g'^ + B.g-""*) 

(77 (B+g-'"" - B^g'"") -7/t (A+g-'^'’ - A_g'*^“) = 

- ^ (B+g-'^ + B-g'"'') 


Omiting the algebraic details, we obtain the relative amplitudes 


irik [T] + k + iCp 12) g-'(^+'?)« 

aI “ 2d 

B_ k-ri[ri-k-iCpl2)e'^’'"'^'>‘‘ 

aI “ 2 B 




irjk 


rf-[k + iCpl2^ 


sin(477a)g 


r-2ika 


2d 

D g-2ika 


where 

and 


A± d 

d = cos(2ria) - sm{2ria)f{k) 


m = i 


k^ + CllA + if 


2T]k 


The equation (37) given us the transmission coefficient 

T = 


1 + sin^(277fl) [|/(A:)|^ - 1 j - 2 sin(4/7fl)Re [/(A:)] 


(31) 

(32) 

(33) 


(34) 

(35) 

(36) 

(37) 

(38) 

(39) 

(40) 


The scattering process is only possible with localized energies in the range IBI > m and then ^ 6 K. We see that 
both the effective potential as the transmission coefficient depends on the energy and then the transmission coefficient 
is not symmetric with respect to E. We can see that the transmission increases with energy, namely, IBI ^ cm T ^ 1. 
The transmission coefficient provides us the necessary condition to obtain the resonant states 


where 77 = 0,1,2,3.. 
2a 


Therefore, the energies of resonant states can be writen as 




and in the limit —> cm, we obtain 


£ = +- 


Nk 

2a 


(41) 


(42) 


(43) 
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For this case r = 0 the energy is fixed {E - C^/Cs - m) and we can see that the pseudoscalar potential influence to 
transmission coefficient as showed in figure 2. From figure 2, we note an oscillatory behavior and no full reflection, as 
expected. As T does not depend on the sing of Cp the transmission coefficient have the same values from both barrier 
and well pseudoscalar potentials. 



c 

p 


Figure 2: Transmission coefficient for r = 0, Ci = 2, m = a = 1. 


Some profiles for transmission coefficient are show in the figure 3 for differents values of Cp and C^. We can 
see that there is not full reflection for all cases. Further, from figure 3 we can note that curves of the transmission 
coefficient for negative energy are closer between them in comparasion to curves for positive energy. As expected, for 
very high energy, T ^ \ and also observe that the pseudoscalar potential maintains the oscillatory behavior of T. 

Also, we can observe that the energies from resonant states is much similar to bound state energies for an infinity 
double-step potential 


ypxi^) = 


—oo 



-t-oo 


X < —a 
X < \a\ 
X > a 


(44) 


in the limit A ^ oo 



(45) 


The similarity is because the bound state energies corresponds approximately to real part of the resonant energies 
obtained from the poles of the transmission amplitude for a square well [31]. A study on correspondence between 
behavior of T and the bound-state energies for nonrelativistic square wells and barriers was done by Maheswari and 
collaborators [32]. But the authors make some mistakes, corrected by Ahmed [33], which also discusses some criteria 
to find oscillatory T for a potential class. 


3.2. Bound States 

The bound-states solutions (IFI < m) can be obtained from (21) using the prescription k - +i\k\ and A+ = = 0 

or k - —i \k\ and - D+ - 0, i.e., the bound-states solutions correspond to the poles of the transmission amplitude. 
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Figure 3: Transmission coefficient for m = a = 1; Cp = 0, Cs = 1 (solid line); = 1, Cs = 0 (solid-dot line); Cp = \, C-z = 0.5 (long dashed 
line); Cp = I, Cz = -0.5 (pointed line). 


The conditions (28) and (29) providing the following quantization condition 


2Ti\k\ 

rj^ + Cl/A-\k\^ 


- tan( 2 ? 7 fl) 


(46) 


We note that the above condition does not depend on the sign of Cp (barrier or square well pseudoscalar potential), 
and hence the energies does not depend on the localization of the delta function. Using the trigonometric relation 
cot( 0 ) = - tan( 0 / 2 ) + ^Jl + coC(O), the quantization condition can be rewrite as 


tan(; 7 «) = + ^|l + lf(k)f - f(k) (47) 

We can see that in the absence of pseudoscalar potential (Cp = 0), the quantization condition is the same obtained 
for spin-0 bosons [34]. This equivalence confirm the results obtained by P. Alberto and collaborators [35], which 
show that the spin and pseudospin symmetries in Dirac equation produce an equivalent energy espectra for relativistic 
spin- 1/2 and spin -0 particles in the presence of vector and scalar potentials. 


3.2.1. Case r > 0. 

For r = 0 we have 77 = +/ \k\, therefore (46) can be writen as 

2\k\^ 

-^ = tanh(2 \k\ a). (48) 

Cl/4-2\kf 

The above condition has just one solution in \k\ - 0, therefore we do not have bound-states solutions for r - 0. 
The case r > 0 has a repulsive barrier between the two delta functions, and does not contain bound-states solutions 
too. The equation (20) and the condition lisl < m allows us to conclude that only have bound states for Cj; > 0. 


3.2.2. Case r < 0. 

For r < 0 we have the condition 0 < Cp < 2mC^ and the quantization conditions provides the figures 4 and 5. 

The effective potential given by equation (19) has not defined parity, and therefore do not know the parity of 
the solutions. Obviously, the case Cp - Q show energy levels with parity defined, even (odd) solutions to negative 
(positive) sign in (47), the same behavior is found in reference [34]. The ground state energy (Eg ^), in the figure 4 
(5), is given by Eg,s - in - C^, {Eg.s - m - CpSgn{Cp)). Therefore we obtain the constrain Cj; -H CpSgn{Cp) - 2m for 
the low-energy state. As we already know, there are always bound-states solutions (at least one) if the intensity of the 
pseudoscalar barrier or well not exceeds the critical value = sj2mC-^, for 0 . 
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Figure 4: Bound-states espectrum for Cp = 0, m = a = 1. 



Figure 5: Bound-states espectrum for C-z = 2, m = a = 1. 


4. Conclusions 

The states of fermions in the framework of mixed vector-scalar-pseudoscalar square potentials was investigated. 
The condition for spin symmetry (A = 0) enables to decouple the Dirac equation in an effective Schrodinger equation 
with a square potential with repulsive and attractive delta-functions situated at the borders for the upper component; 
and the lower component was expressed in term of the upper component in a simple way. An oscillatory transmission 
coefficient and resonant states energy were obtained. We showed that existence of bound-state solutions are condi¬ 
tioned by the intensity of the pseudoscalar potential, which posses a critical value = yfUnC^ for Cj; 0. In the 
absence of pseudoscalar potential, we obtain the same spectrum for spinless particles [34], confirming the predictions 
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of Ref. [35], 

This work can illustrate some general conclusions drawn in previous works about spin and pseudospin symmetries, 
we can obtain the solutions for S = 0 from the A = 0 case, using the chiral transformation (changing the signs of m 
and Cp in the relevant expressions). Finally, it is well known that square potentials, wells and barriers are of certain 
interest in solid state physics, therefore our results could be applied to refine one-dimensional potential models caused 
by ions in a periodic crystal lattice, as the Kronig-Penney model [36]. Other possible application of our results could 
be in the neutron scattering on nucleus, where bound-state informations are extract for many isotopes in well [37]. 
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